function [m b] = minimos_cuadrados(puntosX, puntosY)

    N = length(puntosX);
    sumx = 0;
    sumy = 0;
    sumx2 = 0;
    sumxy = 0;

    for i = [1:N]
        sumx = sumx + puntosX(i);
        sumy = sumy + puntosY(i);
        sumx2 = sumx2 + (puntosX(i)^2);
        sumxy = sumxy + (puntosX(i)*puntosY(i));
    endfor

    den = N*sumx2 - (sumx^2);
    num1 = N*sumxy - sumx*sumy;
    num2 = sumx2*sumy - sumx*sumxy;
    m = num1/den;
    b = num2/den;
endfunction
